Analysis and Numerics of the Chemical Master Equation
نویسنده
چکیده
It is well known that many realistic mathematical models of biological and chemical systems, such as enzyme cascades and gene regulatory networks, need to include stochasticity. These systems can be described as Markov processes and are modelled using the Chemical Master Equation (CME). The CME is a differentialdifference equation (continuous in time and discrete in the state space) for the probability of a certain state at a given time. The state space is the population count of species in the system. A successful method for computing the CME is the Finite State Projection Method (FSP). The purpose of this literature is to provide methods to help enhance the computation speed of the CME. We introduce an extension to the FSP method called the Optimal Finite State Projection method (OFSP). The OFSP method keeps the support of the approximation close to the smallest theoretical size, which in turn reduces the computation complexity and increases speed-up. We then introduce the Parallel Finite State Projection method (PFSP), a method to distribute the computation of the CME over multiple cores, to allow the computation of systems with a large CME support. Finally, a method for estimating the support a priori is introduced, called the Gated One Reaction Domain Expansion (GORDE). GORDE is the first domain selection method in the CME literature which can guarantee that the support proposed by the method will give the desired FSP approximation error. To prove the accuracy and convergence of these three methods, we explore non-linear approximation theory and the theory of Chemical Master Equations via reaction counts. Using these tools, the proofs of the accuracy and convergence of the three methods are given. Some numerical implementations of the three methods are given to demonstrate experimental speed-up in computing an approximation of the CME.
منابع مشابه
Exact results in the large system size limit for the dynamics of the chemical master equation, a one dimensional chain of equations.
We apply the Hamilton-Jacobi equation (HJE) formalism to solve the dynamics of the chemical master equation (CME). We found exact analytical expressions (in large system-size limit) for the probability distribution, including explicit expression for the dynamics of variance of distribution. We also give the solution for some simple cases of the model with time-dependent rates. We derived the re...
متن کاملOn Reduced Models for the Chemical Master Equation
Abstract. The chemical master equation plays a fundamental role for the understanding of gene regulatory networks and other discrete stochastic reaction systems. Solving this equation numerically, however, is usually extremely expensive or even impossible due to the huge size of the state space. Thus, the chemical master equationmust often be replaced by a reducedmodel which operates with a con...
متن کاملApproximate Exponential Algorithms to Solve the Chemical Master Equation
This paper discusses new simulation algorithms for stochastic chemical kinetics that exploit the linearity of the chemical master equation and its matrix exponential exact solution. These algorithms make use of various approximations of the matrix exponential to evolve probability densities in time. A sampling of the approximate solutions of the chemical master equation is used to derive accele...
متن کاملAnalysis of Goldbeter-Koshland Switch Using the Chemical Master Equation
The chemical master equation is considered as an accurate description of general chemical systems, and especially so for modeling gene regulatory networks and protein-protein interaction networks. This paper examines some “off-the-shelf” numerical approaches to solve the chemical master equation for the Goldbeter-Koshland switch on both sequential and parallel machines. The numerical results sh...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
دوره شماره
صفحات -
تاریخ انتشار 2013